Nonlinear stationary solutions of the Wigner and Wigner-Poisson equations 
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Exact nonlinear stationary solutions of the one-dimensional Wigner and Wigner-Poisson equations 
in the terms of the Wigner functions that depend not only on the energy but also on position are 
presented. In this way, the Bernstein-Greene-Kruskal modes of the classical plasma are adapted 
for the quantum formalism in the phase space. The solutions are constructed for the case of a 
quartic oscillator potential, as well as for the self-consistent Wigner-Poisson case. Conditions for 
well-behaved physically meaningful equilibrium Wigner functions are discussed. 
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I. INTRODUCTION 



It is well-known that the solutions for the Vlasov equation are arbitrary functions of the invariants (constants of 
motion) of the system. In the stationary case, this property allows us to construct the so-called Bernstein-Greene- 
^O Kruskal equilibria [ll] for the Vlasov-Poisson system for the classical plasma. On the other hand, it is much more 
^~~^ difficult to derive exact solutions for the corresponding quantum model, namely the Wigner equation or, in the 
self-consistent case, the Wigner-Poisson system. Indeed, even if the classical problem is integrable, the stationary 
r^ ■ Wigner function is not, in general, a function of the classical invariants. This follows since the Wigner time-evolution 
Qh' equation does not preserve the classical constants of motion. Therefore, to date, there is a lack of exact solutions for 
4!) , the Wigner and Wigner-Poisson equations. The better results in this regard are approximate solutions [4|-[j| that 
^ ' have been obtained as the first order quantum correction to the Vlasov-Poisson equilibria. The quantum Bernstein- 
Greene-Kruskal modes have been defined as the solutions of the Wigner-Poisson system under the periodic boundary 
conditions [a], but the explicit construction of these solutions is still a challenge. Notice that the original article by 
Wigner [y|, where he introduces his celebrated function, considers the first-order quantum correction, to a Maxwell- 
Boltzmann thermodynamic equilibrium. In addition, quantum like corrections were proposed for charged-particle 
\^ ' beam transport 0]. 

The purpose of the present paper is to develop explicit exact solutions of the Wigner equation, in a particular form 
that is described in Section II. Both the external and self-consistent potential cases are treated. In a sense, we develop 
Qv a quantum analogue of the Bernstein-Greene-Kruskal equilibria in the phase space. Of course, exact stationary Wigner 

(T^ ■ functions can be found after applying the Wigner transform to the previously known energy eigenstates, if available. 

, I I However, the question is to derive solutions for the Wigner equation considered in itself, without reference to the 

f^ • Schrodinger equation or the nature of the quantum statistical mixture of the system. Here we propose a special 

QQ functional Ansatz, as a function of the energy and the position, for the Wigner function in a conservative system. The 

(T^ , technique is illustrated for a quartic anharmonic oscillator case and for the self-consistent Wigner-Poisson system. 

L| ■ The derived solutions are the first explicit exact solutions of the Wigner equation for nonlinear systems, constructed 

• i-H , independently of the Schrodinger formalism. In this way, we can have insight on the relation between the classical 

p\^ ' constants of motion and the solutions of the stationary Wigner equation. In this regard, it is already well known 

JH \ [8| that for a quadratic Hamiltonian there is a close relationship between the Ermakov-Lewis invariant, which is the 

■ - - ' basic classical constant of motion of the problem, and the Wigner function. However, for a quadratic Hamiltonian the 

Wigner equation reduces to the classical (Vlasov) equation, with the quantum effect restricted to the initial conditions. 

The work is organized as follows. In Section II, we consider a quartic nonlinear oscillator problem and derive 

quantum phase space structures by using our Ansatz generalizing the Bernstein-Greene-Kruskal solutions for the 

classical system. In Section HI, the approach is adapted to the Wigner-Poisson system, which is related to a quantum 

plasma system. Section IV is reserved to the conclusions. 



* Also at Universidade do Vale do Rio dos Sinos - UNISINOS, Av. Unisinos 950, 93022-000, Sao Leopoldo, RS, Brazil 

t Also at Department of Physics, Umea University, SE-90187, Umea, Sweden; GOLP / Instituto de Plasmas e Fusao Nuclear, Instituto 

Superior Tecnico, Universidade Tecnica de Lisboa, 1049-001 Lisboa, Portugal; SUPA, Department of Physics, University of Strathclyde, 

Glasgow, G40NG, UK; School of Physics, University of Kwazulu-Natal, Durban 4000, South Africa. 



V 



V 





X 



FIG. 1: the potential in Eq. ([6} for /i = 1 and /i = —1, respectively. 



II. EXACT SOLUTIONS OF THE STATIONARY WIGNER EQUATION FOR A NONLINEAR 

OSCILLATOR POTENTIAL 



The Wigner equation [9'| in one spatial dimension with a potential V{x,t) reads 
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V{x+—,i)-Y(x-—,t) 



f{x,v',t), 



where f{x, v, t) is the Wigner pseudo-probability distribution and all quantities have their usual meaning. The formal 
classical limit {h -^ 0) of the Wigner equation is the Vlasov equation, with the Wigner function playing the role of the 
probability distribution function. Even assuming negative values in some regions of phase space, the Wigner function 
can be used to compute the macroscopic quantities like the density and the current, in the same way as a faithful 
distribution function. In addition to Eq. ([T]), a genuine Wigner function should correspond to a positive definite 
density matrix. Therefore, / = f{x,v,t) must satisfy J^] at least the following necessary conditions. 



dxdvf — 1 



dvf > 0, 



dxf > 0, 



dxdvf < 
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(5) 



Equation ^ is just a normalization condition, while Eqs. ([3]) and ^ arise because the spatial and velocity densities 
should be non-negative everywhere. Finally, Eq. ([5]) is needed to avoid violation of the uncertainty principle, ruling 
out a too spiky function f{x,v,t). 

In the stationary case, one would be tempted to use the conservation of the energy to construct exact solutions for 
Eg. ID). However, except in the somewhat restricted case where the Wigner function is a linear function of the energy 
[llj, the Wigner equation is not satisfied by functions of the energy alone. To illustrate our technique, we consider 
the anharmonic potential 
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where /z = ±1 is a numerical parameter and u and k are parameters with dimension of angular frequency and 
wavenumber, respectively. Plots of the potential are shown in Figure 1, showing a single nonlinear potential well 
or a double-well potential, according to the values of /i. The coefficients were chosen to match the expansion of a 
pendulum-like potential ^m{uj'^ / k'^)[l — cos(fca;)]. A similar form was considered in connection with the numerical 
simulation of quantum echoes described by the Wigner equation [12| . 



For the chosen potential, in the stationary case {df/dt = 0), the Wigner equation ([T]) is 



df , 2 
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If there was no nonhncarity, Eq. ([7]) would be the stationary Vlasov equation for the simple harmonic oscillator whose 
solutions are arbitrary functions of the energy. However, in the quantum and nonlinear case, the higher-order velocity 
derivative term prevents the existence of a solution as a function of the energy only. Notice that the term in the 
right-hand side of Eq. ([7]) can assume large values: it is not necessarily just a quantum correction. Higher-order (in 
h ) terms would appear for higher anharmonicities. 

With the rescaling q — kx.p — kv/uj. F — ujf /k"^^ Eq. ([7]) transforms into 



where 
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is a non-dimensional parameter measuring the relevance of the quantum effect. At this point, notice that the necessary 
conditions (O^® are rewritten as 



dqdpF — 1 , 
dpF > 0, 
dqF > 0, 
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27r\/6r 

It is more convenient to search for solutions of Eq. fS]) in the form F = F{H, q), where 

^ = y+ny-24 



is the energy. Then from Eq. ([5]) we obtain 
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The usefulness of the energy and position variables in the treatment of the Wigner equation has been recognized in 
the semi-classic case (3|, but we feel that this method can be pursued in more depth for a fully quantum system. 
When there are no quantum effects, the right-hand side of Eq. ([15]) is zero and the solution is simply F — F{H), 
with an arbitrary functional dependence, in the same spirit of the Bernstein-Greene-Kruskal solutions for the classical 
plasma. The difficulties with Eq. (|15p are in the entangled character of the quantum term, which prevents the use of 
the separation of variables technique, for instance. It is also apparent that if dF/dq = 0, then any linear function of 
the energy would be a solution. 

Our proposal is to consider exact solutions according to 



F=[A{q)H + B{q)]e 



C{q)H 



(16) 



where A, B and C are functions to be determined, depending on position only. For A = 0, the Wigner function 
automatically has a Gaussian shape in velocity space. Otherwise, for A ^ 0, one can have a two-stream type Wigner 
function, double-humped in velocity space, at least for specific parameters. Notice that the proposed solution does 
not distinguish between the trapped and untrapped particles. This is another different feature in comparison to the 
Bernstein-Greene-Kruskal method. 



After inserting Eq. (fT6|) into Eq. (|15p . the exponential factorizes and we derive a second-degree polynomial of H, 
which must be identically zero. Setting the coefficients of the different powers of the energy to zero, we obtain 



A' + BC' = liVq 
B' = fiTq 
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where the prime denotes derivative with respect to q. 

For i? = 0, it can be verified that the above system do not admit any solutions. For B y^ 0, the result is 
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where Aq,Bq and Cq ^ are integration constants. If Aq = 0, then the Wigncr function is certainly Gaussian in 
velocity space. 

In the double- well potential case, when /i = — 1, the exact solution is real and bounded provided Co > 0. However, 
for /u = 1, it can be shown that the exact solution becomes singular when g2 _^ Co/(2r), unless < Cq < 24r. 
Actually, when n — 1 the Wigner function in Eq. (fT6|) is acceptable only for q^ < Co/(2r). Figure 2 displays a 
typical plot in phase space, for F = 1,/i = l,^o = O-.Bq ^ 0.42, Co = 16. For simplicity, we assume that there 
are no scattering states, so that F = for q^ > Co/(2r) — 8. The Wigner function has a Gaussian shape in the 
momentum space and an abrupt localization in configuration space, as can be seen more clearly in Figure 3 for the 
same parameters and for different rescaled velocities p = and p = 2. These strongly localized structures do not 
have classical counterpart. Observe that the necessary conditions in Eqs. pHll - p^ are also satisfied. In particular, 
JdqdpF^ = 0.03 < (27rV6r)-i = 0.07. 

When /i = —1, the Wigner function is regular provided Co > 0. Then, we can have a more rich variety of behaviors. 
For instance, for fi — —1, Aq > 0, Bq — — 2.5Ao, Co = 1, there is a two-humped distribution in velocity space, provided 
F > 0.37. Figure 4 shows the Wigner function in momentum space, for fixed position q — 0, with the same parameters 
and several values of F. One sees that for increasing quantum effects there is a progressive depth of the exact solution, 
which eventually becomes negative for F > 0.63. Notice that, in principle, there is no limiting value of F, since the 
solution is non perturbative. However, for a physically meaningful positive particle spatial density 
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there are extra restrictions, even recognizing that in general the Wigner function is not a positive definite quantity. 
For instance, when /i = —1,^0 > 0;^o ~ ^2.574o,Co ~ 1, n{q) > for all q imposes that F < 0.3. For these 
parameters, the other necessary conditions in Eqs. pO|) . (fT2|) and (|13p are also fulfilled, as can be checked after 
numerically performing the spatial integrations. Actually, a detailed calculation shows that when ^ = — 1, ^o > and 
Bq < 0, for a positive definite particle density the exact Wigner function cannot be a two-stream type distribution. 
On the other hand, for /i = — 1 and non negative ^O: Bq (in which case the distribution function is not two- humped 
in momentum space), in principle one can have arbitrarily large values of F. 

The opposite behavior arises when fi = —1 and Aq < 0, i?o > 0. In this case, it can be proved that a positive 
definite particle density implies a two-humped Wigner function in momentum space. For instance, at the origin g = 0, 
the condition for n(0) > reads 



UBqCq^^T^ > \Aq\ (c2 + 6FCo + 24F^ 



(22) 
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FIG. 2: the exact Wigner function in Eq. ^ for T = 1, /i = 1, Ao = 0, Bo = 0.42, Co = 16. 




FIG. 3; spatial localization of the exact Wigner function in Eq. (|16p for F = 1, /x = 1, Ao = 0, -Bo = 0.42, Co = 16 and p = 
(solid line) and p = 2 (dashed line). 



and this inequality can be shown to imply that the Wigner function admits two minima in velocity space, as illustrated 
in Figure 5. Notice the opposite concavity in comparison to the /i = —l,Ao > 0,Bo < case. In the same figure, 
we show the particle density, which has a deep centered at the origin but is positive definite, even with the Wigner 
function admitting negative values. 

Now, it can happen that the quantum effect comes in favor of a physically meaningful solution, since n{q) > 
for all q provided F is sufficiently large. For instance, from Eq. (P^ it follows that n{0) > if BqCq > 2|Ao| and 
F is large enough. In particular, for fi — ~l,Ao < 0,Bo — —2.5Ao and Co — 1, one needs F > 1.14. Figure 6 
shows n{q) for these parameters and several values of F. One sees that the two maxima of the particle density tend 
to approximate each other as the quantum parameter increases. Also, the minimum and maximum values of n{q) 
increases as quantum effects becomes larger. However, F can't be arbitrarily large, since the solution becomes too 
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FIG. 4: the solution in Eq. (|16p in momentum space for q = 0, fi = —1, ^o ~ 0.52, Bo = —1.3 and Co ~ 1. The quantum 
parameter is F = 0.2 (left, top), F = 0.4 (right, top), F = 0.6 (left, bottom) and F = 0.8 (right, bottom). In the graphs, only 
for F = 0.2 the solution is a physically acceptable Wigner function. 



localized in phase space for increasing F. For instance, when F = 5,/j, = — 1, Aq = —0.11, Bo = 0.28 (figure 6, right, 
bottom), we have J dqdpF^ = 0.032 > (27r-\/6r)~^ = 0.029, violating Eq. (fT3| . In conclusion, in each specific case 
the necessary conditions (fTU)) - p^ must be checked. 



III. APPLICATION TO THE WIGNER-POISSON SYSTEM 



Now consider the self-consistent case where the potential is V 
the electrostatic potential satisfying the Poisson equation 



-ecp, where e is the elementary charge and (j) is 






= "(</') = fdv, 



(23) 



where eo is the permittivity constant and ng a background ion number density. Coupled to the Wigner equation, 
Eq. p3p compose the Wigner-Poisson system [1J|. For the latter, it is convenient to search for solutions in the form 
/ = f{H, n), where H — mv^ jl — ecf) is the energy and n is the electron number density, as defined in Eq. (|23p . Then, 
writing the stationary Wigner equation, retaining only the first-order quantum correction, 



we obtain, after using the Poisson equation 
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FIG. 5: on the top: the Wigner function in Eq. ([161 
bottom: the corresponding particle number density. 
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FIG. 6: particle number density in Eq. 121} for ^ = -1, So = -2.5Ao, Co = 1 and for F = 0.5, Aq = -0.48 (left, top), T = 1.14, 
Ao = -0.13 (right, top), T = 2, Aq = -0.108 (left, bottom) and T = 5, Aq = -0.111 (right, bottom). For each F, a different Aq 
must be chosen, to comply with a normalized Wigner function. The derived solution complies with the conditions P0|) - (I13|) . 
with the exception of the case F — 5, Ao = —0.111. 



where Up = (noe^/meo)"^/^ is the electron plasma frequency. In Eq. (^5]) . the scalar potential is interpreted as a 
function of the electron number density. This is locally possible using n = n{(f)), according to the implicit function 
theorem, provided n is not identically a constant. Now the Ansatz 

f = [A{n)H + B(n)]e^(")^ (26) 

produces the following system for the functions A, B and C depending on the electron number density only, 
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where the prime denotes derivative with respect to n. The resulting system can be solved in all generality, but here 
we are content with the A = case, which is more amenable to detailed calculations. When A = 0, the solution for 
the system ([27|) is given recursively as 
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(1 - n^ujlC^n/{6no)y/^ ' 



Bo I 1 



6no 



-3/4 



exp 



12no 



dn4){n)C^{n) 



(28) 
(29) 




FIG. 7: the electrostatic potential as a function of the electron number density, for e = m = Ti 
(upper curve, dot-dashed), Co = 2 (mid curve, dashed) and Co = 3 (lower curve, solid line). 
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Bo = 1 and Co 



where Bq and Co are integration constants, with Co > for an intcgrablc Wigncr function. 
Using Eq. (|^5|) . we obtain 



n{(j)) = B 
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After rearranging, Eq. ([M)) then gives the foUowing integral equation for B — B{n) 
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exp 
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where C — C{n) is given by Eq. 

In the absence of the quantum effect, from Eq. (|29p . one would directly have B = Bq. Since in the self-consistent 
case we are considering, for simplicity, only the first-order quantum correction, it is then reasonable to put B = Bq 
in the integral at the right-hand side of Eq. pip . After performing the integral, which is a complicated expression 
involving the dilogarithm function Li2(2:) = L dshi{s)/{l — s), one derives a function B = B{n) and hence the 
electrostatic potential = 0(n) through Eq. (|30)) . As apparent from Eqs. (|28|) . (|29|l and ([3T]) . one can identify 
the non-dimensional quantum parameter Ti uj'iC'^/Q, with the interpretation of the inverse square of a temperature. 
Figure 7 displays typical plots of (j){n), with units so that e = m = h = uq — ujp = Bq = I, for different values of 
Co. The graphs are similar to that from a electron hole Maxwell-Boltzmann equilibrium (0 ^ Inn), however with 
modifications and a cutoff when n -^ 6no / {h^ uj^Cq) . For increasing quantum effects (larger Cq, or a smaller effective 
temperature), the cutoff becomes closer to the origin. The explicit form in terms of position ((/> — 0(x)) could also be 
found, through the Poisson equation. 

In the Wigner-Poisson case, the calculations becomes too involved if higher-order quantum corrections are added in 
the right-hand side of Eq. ([M]) . However, the stationary solution found are non perturbative, in the sense that they 
do not have, as a starting point, a prescribed zeroth-order Wigner function, like a Maxwell-Boltzmann or Fermi-Dirac 
equilibrium, as in Refs. [4]-[6[. In other words, except for the substitution B ^ Bq in the right-hand side of Eq. (|3ip . 
what have been found here is an exact solution for an approximate model, namely the Wigner-Poisson system up to 
the first-order quantum term. 



10 

IV. CONCLUSION 

In this paper, we have presented exphcit nonhncar solutions for the stationary Wigner and Wigncr-Poisson equa- 
tions. The solutions are inspired by the Bernstein-Greene-Kruskal modes of the classical Vlasov-Poisson plasma. 
However, unlike for the Bernstein-Greene-Kruskal modes, the equilibrium distribution function here depends not only 
on the energy, but also on the position. In addition, there is no need for the eigen-functions or the characterization of 
the quantum statistical mixture of the system, with no reference to the Schrodinger equation. Our new solutions are 
in the form shown in Eqs. (|16p or (|26p . which disentangle the equilibrium Wigner equation, thanks to the particular 
exponential dependence on the energy. It would be relevant to obtain another classes of explicit solutions, in order to 
derive a better understanding of the relation between the classical constants of motion and the equilibrium states of 
the Wigner and Wigner-Poisson equations. Finally, there are no known necessary and sufficient criteria for a faithful 
Wigner function in the phase space. The solutions derived have been checked against the necessary conditions Eqs. 
(HI)-®, for specific parameters only. 
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